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Abstract 

We numerically analyze the random matrix ensembles of real-symmetric matrices with col¬ 
umn/row constraints for many system conditions e.g. disorder type, matrix-size and basis- 
connectivity. The results reveal a rich behavior hidden beneath the spectral statistics and also 
confirm our analytical predictions, presented in part I of this paper, about the analogy of their 
spectral fluctuations with those of a critical Brownian ensemble which appears between Poisson 
and Gaussian orthogonal ensemble. 
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I. INTRODUCTION 


Random matrices with column/row constraints appear as the mathematical represen¬ 
tations of complex systems in diverse areas e.g bosonic Hamiltonians such as phonons, 
and spin-waves in Heisenberg and XY ferromagnets, antiferromagnets, and spin-glasses, eu¬ 
clidean random matrices, random reactance networks, hnancial systems and Internet related 
Google matrix etc (see section H and references in m)- It is therefore relevant to seek the 
information about their statistical behavior; this not only helps to probe the related com¬ 
plex system but also in revealing the hidden connections among seemingly different complex 
systems. 

The standard route for the statistical analysis of a random matrix ensemble is based 
on its eigenvalues/ eigenfunctions fluctuations. This in turn requires a knowledge of the 
joint probability distribution function (JPDF) of the eigenvalues and eigenfunctions. An 
integration of JPDF over undesired variables then, in principle, leads to the fluctuation 
measures for the remaining variables; for example, the spectral fluctuation measures can be 
obtained by an integration over all eigenfunctions. As discussed in [T], the JPDF for the 
column/row constrained ensembles can be derived from the ensemble density but its further 
integration over eigenfunctions is not an easy task. This motivates us to seek, an alternative 
route i.e the numerical investigation of the ensembles; our focus in this paper is on following 
objectives: 

(i) to understand the effect of randomness (both type and strength) and basis-connectivity 
on the fluctuation measures. 

(ii) to analyze the size-dependence of the fluctuation measures, 

(iii) search for the existence of any universal features in the statistical behavior. 

The motivation for the first objective comes from previous studies of the systems with 
Goldstone symmetry; the column constraints in the matrix representation of their operators 
originate from the invariance under a uniform shift pQ. In absence of disorder, the low fre¬ 
quency excitation in these systems have long wavelengths and are described by the equations 
of motion based on macroscopic properties. Disorder causes local fluctuations in the prop¬ 
erties but excitations are believed to couple only to their averages over a volume with linear 
dimensions set by a wavelength. Gonsequently, higher dimensions, where averaging is more 
effective, are expected to be disorder-insensitive HD]. But almost all these results are based 
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on average behavior of the spectral density and inverse participation ratio (IPR). As the 
statistical properties can in principle be expressed in terms of the eigenvalues/eigenfunctions 
fluctuations of the related operators, it is natural to query whether disorder has any, and 
if so, what impact on the fluctuations of the excitations, especially in higher dimensions. 
As mentioned in section 2 of [1], the spectral statistics of the low lying bosonic excitations 
above a ground state Eq is analogous to that of an ensemble of column constrained matrices 
with column constant ai = Eq (for Z = 1, 2, ..A^). The latter being statistically similar to the 
case a; = 0 (except for a shift in the spectrum), the spectral information about Goldstone 
modes can be obtained by an analysis of the = 0 ensemble. Keeping this in view, the 
present analysis is conhned to the zero-column constraint ensembles. 

The modelling of complex systems by random matrix ensembles usually requires consid¬ 
eration of the inhnite size limit. The numerical analysis however involves only hnite size 
matrices. To translate the numerical information to real systems, it is necessary to analyze 
hnite size effects which forms the basis of our 2nd objective. 

The necessity for the third objective originates from our recent study |T] of the JPDF 
of the column constrained ensembles (CCE) for two different disorder types and basis- 
connectvities. Our results indicated a mathematical analogy of their JPDF with that of 
a Brownian ensemble (BE) appearing between Poisson and Gaussian orthogonal ensemble 
(GOE); this in turn implies a GGE-BE analogy for the local spectral huctuation measures. 
The implication is seemingly in contrast to a previous study [1] which predicts the analogy 
of the two point level-density correlation of GGE to that of a GOE. But, as explained in pQ, 
the BE-statistics in an energy range near e = 0 is indeed of GOE type (see the discussion 
in 2nd paragraph after eq.(73) in [1] ) and, therefore, our prediction agrees with |1] near 
e = 0. For other energy ranges, a direct numerical comparison of their fluctuation measures 
is needed to reconhrm the GGE-BE analogy. (Although a comparison of GGE numerics 
with analytical BE-results would have been more desirable but, due to technical complexity 
related to the orthogonal-space integration, available BE-results are mostly asymptotic). 

To fulhll the above objectives, here we numerically study GGEs for many system con¬ 
ditions. More specihcally, we consider the Gaussian and bimodal distributions of the off- 
diagonal matrix elements for many disorder strengths and system sizes; the diagonals are 
obtained by invoking the column sum rule. The dimensionality in the matrices is taken into 
account by varying basis-connectivity i.e by choosing the matrix elements forming nearest- 
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neighbor pairs non-zero and all others zero (still subjected to column sum rule). 

The paper is organized as follows. Section 2 presents the ensemble densities for the cases 
considered in our numerical analysis; these ensembles are special cases of those considered in 


section 4 of [T] . An analysis of the spectral fluctuation measures requires a specihc rescaling 
of the spectrum which in turn depends on whether the spectrum is ergodic or not; this is 
discussed in section 2. A alongwith a numerical study of the ergodicity for CCEs, BEs as well 
as GOE; the latter is included to illustrate the deviation of CCEs from GOE. Sections 2.B and 


2.C present the numerical results for the average behavior and fluctuations, respectively, of 
CCEs with Gaussian and bimodal distributed off-diagonals, both cases considered for inhnite 
as well as hnite basis-connectivity and many matrix sizes; the objective of these sections is 


not only to verify the theoretical insight derived in [T] but also to provide information about 


mathematically intractable fluctuation measures. The results of a numerical comparison of 
the CCE-BE fluctuations are given in section 3 which reconhrm our theoretical predictions 
of pp. Section 4 concludes this study with a brief reviews of our main results and open 
questions. 

II. COLUMN CONSTRAINED ENSEMBLES: FOUR TYPES 

To investigate the sensitivity of the statistical fluctuations to the nature of randomness 
and basis-connectivity, we consider following CCE ensembles of real-symmetric matrices H 
with p{H) as their probability densities: 

(a) independent off-diagonals with Gaussian distribution (infinite range or d = oo CCGE 
case) 


This case corresponds to case II of pQ with a = 0. 



( 1 ) 


(b) Independent, nearest-neighbor Gaussian hopping in d-dimension (finite d GGGE case) 
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with z as the number of nearest neighbors of a site. The symbol to a sum over 

all nearest-neighbors of a given site k (excluding I = k site). This case corresponds to case 
III of [I] with bo = bi = 0, a = 0. 


(c) Independent off-diagonals with himodal distribution (inhnite range or d = oo CCBE 


case 


N \ ^ 

p(n)^[ n wwi-1 )+hhm +1)]] 


kl 


( 3 ) 


1=1 


This case corresponds to case IV of [T] with a = 1, a = 0. 

(d) Independent, nearest-neighbor himodal hopping in d-dimension (hnite-d CCBE case) 


N 


p(H)=N I n 

.k,l;lGz{k) 


S{IIki - qa) 

.g=±i 



N 


N 


k,l-,l^z(k) 


kl 


1=1 


( 4 ) 


This case corresponds to case V of pQ with a = l, 6 i = 0,a = 0. 

For numerical computation of the eigenvalues and eigenfunctions of each of the ensembles 
above, we exactly diagonalize thousands of matrices of sizes N = 512,1000, 2197,4913 using 
LAPACK subroutines; by choosing an appropriate number of matrices, the statistical data 
size for each N is ensured to be same. Further, as clear from eqs.( Tppp ), the spectral 
statistics for these cases is independent of the off-diagonal variance (a rescaling of the matrix 
elements by the variance leaves the statistics unaffected); this is also confirmed by our 
numerics (shown here for the level-density case in hgure 1). For fluctuations studies, the 
variance of non-zero off-diagonals in eqs.(T]| 2 ) is therefore kept fixed at 7 = 1 / 2 . 

Prior to calculation of the fluctuation measures, it is necessary to determine the ergodic 
as well as stationary aspects of the spectrum. 


A. Ergodicity and stationarity of the spectrum 

The spectral density Pe(e) at a spectral scale e is dehned as Pe(e) = ~ ^k)- 

Under assumptions that pe consists of two spectral scales, a slowly varying component psm, 
often referred as the secular behavior, and a rapidly varying component, pjiuc, one can write 
Pe = psm + Pfiuc (also referred as the local random excitations of p^m) [2113]. A spectral 
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averaging over a scale larger than that of the fluctuations gives Pfiud^) = 0 and 

Psm = ^ Ie-Ae /2 Pei^)■ Previous studies of a wide range of complex systems indicate psm 
to be more system specihc as compared to the fluctuations; the latter reveal a great deal 
of universal features and often depend only on the global constraints of the system. The 
information about local fluctuations can therefore be better revealed by an ’’unfolding” the 
spectrum i.e removing its system-dependent part psm by a rescaling = J^'^Psm{x) dx of 
the eigenvalues which results in uniform local mean spectral density Eli- 

Experimental studies often analyze the behavior of a single spectrum which leaves the 
unfolding by psm the only option. But theoretical modelling of a complex system by a random 
matrix ensemble is based on the assumptions of ergodicity of global and local properties. As 
discussed in section 2.2.3 of [6|, the ergodicity of a global property such as the level density 
implies Psm{^) = -Ri(e), with Ri{e) as the ensemble averaged level density, and one can use 
i?i(e) as a substitute for Psm{^) for various analytical purposes. It is therefore necessary to 
check the ergodicity of the level density which requires mm 

{Psm{ef)-{psm{,e))^^0, {psm{e)) = Ri{e) (5) 

with (.) implying an ensemble average for a hxed e. A comparison of the ensemble and the 
spectral averaging for the spectral density of CCGE, BE as well as GOE is shown in hgure 
1; here each case is considered for an ensemble of 5 x 10^ matrices of size N = 1000. As seen 
from the hgure 1(d), Psm{e) even for a single GOE matrix almost coincides with {psm{e)) as 
well as Ri{e) for all e, reconhrming its ergodicity. For GGGE and BE, the e-dependence of 
the psm huctuates (weakly) from one matrix to the other but the ensemble averaged Psm{^) 
does agree well with Ri{e). Although the hgure 1 displays the result only for one system size, 
we hnd similar analogy for higher N values too, with deviation from one matrix realization 
to other diminishing with higher N. This indicates an approach to ergodicity for the level 
density in large N limit. 

An ergodic level density however does not by itself implies the ergodicity of the local 
density huctuations. As discussed in section 2.2.3 of [B], the ergodicity of a local property 
say ”/” can be dehned as follows: let Pe{s) be the distribution function of / in the region 
”e” over the ensemble and Pk{s) be the distribution function of / over a single realization, 
say ”fc” of the ensemble. / is called ergodic if it satishes the following conditions: (i) Pe{s) is 
stationary i.e Pe{s) = Pe'{s) for arbitrary e,e', (ii) for almost all k, Pk{s) is independent of k 
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i.e Pk{s) = P{s), (iii) Pe{s) = P{s). However if condition (i) is not fulfilled, / is then referred 
as locally ergodic. More clearly, an ensemble is locally ergodic in a spectral range say Ae 
around e, if the averages of its local properties over the range for a single matrix (referred 
as the spectral average) are same as the averages over the ensemble at a hxed e (referred as 
the ensemble average). To analyze the ergodicity of local spectral fluctuations, we consider 
a standard measure, often used for this purpose, namely, the width a{k, e) of the distribution 
P{k;s,e) of the eigenvalue spacings Sk between k^^ neighbors (k = 0, l,2..j within a range 
(e — Ae/2, e + /S.e/2) with s = Sk/D and D as the mean spacing [H E]; (note the mean of 
the distribution is k + 1). Figure 2 compares the ensemble averaged ak with the spectral 
average for CCGE, BE and GOE, separately for each case (using eq.(35) and eq.(39) of |6] 
for the spectral and the ensemble average, respectively)). As the hgure indicates, contrary to 
GOE, the local fluctuations in both GGGE as well as BE are non-ergodic and their unfolding 
with i?i(e) for hnite N cases is not appropriate (as it may introduce spurious fluctuations). 
It is therefore necessary to determine the local spectral density psm either by theoretical, 
experimental or computational route. 

As mentioned above, the stationarity is an important aspect which implies an invariance 
of the local fluctuation measures along the spectrum-axis. A non-stationary spectrum can at 
best be locally ergodic and therefore a complete information about the spectrum requires its 
analysis at various spectral ranges. Figure 3 displays the behavior of an ensemble averaged 
ct(/c; e) for GGGE, BE and GOE for many fc-values and for a few spectral ranges taken 
from the region between edge and bulk. Note, as already known [5], GOE is also not 
exactly stationary and the edge statistics shows signihcant deviation from that of the bulk 
(our numerics shown in hgure 3(d) reconhrms this behavior). As clear from the hgure, the 
regions in which translational invariance holds for GGGE and BE are much smaller than in 
GOE. 

B. Spectral density and inverse participation ratio: average behavior 

The average spectral density psm is a standard measure for experimental or computational 
studies of many complex systems (often based on single system analysis); hgures 4,5 depict 
Psm behavior of a column constrained matrix, for two disorder types, namely, Gaussian and 
bimodal, for various dimensions (i.e basis connectivity) and disorder strengths. The ehect of 
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dimensionality on the average behavior is clearly visible from these hgures. As shown for the 
Gaussian case in hgure 4, the dependence of psm on the disorder strength 7 can be scaled out. 
The size N has no effect on ps^ for short-range basis connectivity (nearest-neighbor hopping 
cases with low-dimensionality d < 3), for both type of disorders (ensemble ( [2pl )) . The 
numerics for different N for the long-range basis-connectivity (ensemble 1 ®). displayed in 
hgure 5(c), however indicates following A^-dependence Psm{e) = -^f ^ comparison 

shown in hgure 12 (a) suggests the bulk behavior of the function / as semicircle : Psm(e) = 
-^\/2aN — with a a constant. An important point to note from hgure 5(c) is the same N 
dependence of psm in both edge and bulk; the critical BE with /i oc A^ is expected to model 
the d = 00 CCE-statistics in both energy regimes (see section V.B of P for the details). 
Furthermore the edge-behavior seems to be in agreement with the level-density of Goldstone 
modes in spin-glasses for d > 3: Psmifjj) ^ where u = e — cq and g = 2 for d = 3 and 
q = 1.5 for d = 00. (Based on various studies, the behavior of level density for Goldstone 
modes in spin glasses is believed to be as follows: Psm{<-^) cc for d = 1 , Psm{^) cc oo'^ for 

d = 3 mm and Psm{’^) oc for inhnite-range spin glass [ini[i3]). Here, for clarity of 
presentation, the behavior for d?i(e) and {psm(e)) are not included in hgures 4,5 but, based 
on hgure 1 , their dependence on disorder and system-size is expected to be same as that of 

psm- 

As discussed in P, a typical eigenvector of a real-symmetric column constrained matrix 
displays a new characteristic: it is the lack of single basis-state localization (the minimum 
number of basis-states required for localization is two) which is in constrast with an uncon¬ 
strained real-symmetric matrix. This tendency can be better analyzed by a numerical study 
of the inverse participation ratio (IPR) I 2 , the standard tool to describe the localization 
behavior of an eigenfunction. It is dehned as hiOn) = J2k=i for an eigenfunction On 

with Okn as its components in a At-dimensional basis. The minimum number of basis-states 
available for localization restricts the maximum IPR /2 < | with k = 1 or 2 for a typical 
eigenvector of a matrix with or without column constraints, respectively P . As can be seen 
from hgure 6 , the maximum IPR indeed does not go above 0.5, except in hgure 6 (f) where 
it reaches 0 . 6 ; the latter seems to be an edge ehect. 

The localization length in general huctuates from one eigenfunction to the other. The 
average localization behavior of the eigenfunctions within a given spectral range can be given 
by a spectral averaged I 2 (i.e average of IPR of all eigenfunctions within the range) [ 12 ]. In 



case of an ensemble, it is usefnl to consider the ensemble as well as spectral averaged I 2 , later 
referred simply as ( 12 ). Fignre 6 indicates the existence of a scale-invariance of {I 2 ) with 
respect to off-diagonal variance and therefore degree of disorder for both types. As clear 
from the insets in hgnre 6, the localization tendency for the low lying eigenfnnctions (lower 
energy excitations) for the Ganssian and Bimodal types is jnst the reverse for short-range 
basis connectivity bnt similar for the long range. The hts in the insets for each case snggest 
{l2,sm) ~ for d = 1, {l 2 ,sm) ~ for d > 3 for small-o;. Previons spin glass stndies on 
the Goldstone modes give the avearge (spectral) localization length (^( 0 ;) oc for d = 1 

and C(i^) oc for d > 3 [lU] which implies, nsing l 2 ,sm oc [IB], h oc for d = 1 
and l 2 ,sm oc for d > 3. As the hts in the inset indicate, the GGE-resnlts are in agreement 
with spin glass behavior. 

C. Spectral fluctuations 

Our next step is to probe the local fluctuations around the spectral density for two 
disorder-types and for different connectivity in the lattice. As discussed above, Psm{e) for 
GGE is e-dependent and a comparison of local fluctuations from two different spectral re¬ 
gions with different Psm{e) requires an unfolding of the spectrum. In past, many ways 
for unfolding have been suggested e.g. polynomial unfolding (approximating the empirical 
staircase function by a polynomial) [7|, Gaussian unfolding (Gaussian smoothing of each 
delta peak in pe) EE], Fourier unfolding (removal of short scale wavelengths using Fourier 
transformation) or local unfolding (local mean level density assumed to be constant within 
a narrow spectral-window) [8], detrending the fluctuations [9] etc. Previous studies of many 
systems indicate the statistics in the low-density regions e.g. edge to be sensitive to the 
procedure used but that in the bulk is unaffected. Although Gaussian unfolding is often 
believed to be suitable in the edge region, the results are sensitive to width of the Gaussian 
smoothing [8]. Our numerical experiment with various unfoldings encourages us to apply 
the local unfolding [8|: we first determine the smoothed level density psm for each spectra by 
a histogram technique, and then integrate it numerically to obtain the unfolded eigenvalues 
Tn = Psm de. (Note, this is similar to the Gaussian unfolding except now the smoothing 
of ^-functions is done ove the bins of fixed size). Testing of this unfolding procedure for two 
fluctuation measure of a GOE, namely, the nearest neighbor spacing distribution and the 
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number variance gives results consistent with the GOE-theory, for both bulk as well as in 
the edge regime. Note, due to spurious fluctuations introduced by hnite size-effects, all the 
unfoldings mentioned above are unreliable for large spectral length-scales. More clearly, the 
numerical calculation of the short-range fluctuation measure is robust across different un¬ 
folding procedures but long-range measure are sensitive to the type of unfolding used. Only 
way to distinguish the ’’true” fluctuations (those present in iV —)■ cx) limit) from the spurious 
ones is to analyze the statistics for many system-sizes and study the emergent behavior in 
large size limit [5j. 

The fluctuations being in general sensitive to the spectral-range (due to non-stationarity 
of the spectrum), it is necessary to choose an optimized range AE which is sufficiently large 
for good statistics but keeps mixing of different statistics at minimum. We analyze 2-5% 
of the total eigenvalues taken from a range AE, centered at the energy-scale of interest 
(edge or bulk). This gives approximately 10® eigenvalues and their eigenfunctions for each 
ensemble. (Note due to almost constant level density in the bulk, the statistics is locally 
stationary and one can take levels within larger energy ranges without mixing the statistics. 
A rapid variation of psm in the edge however allows one to consider levels only within very 
small spectral ranges. For edge-bulk comparisons, it is preferable to choose same number of 
levels for both spectral regimes). 

The standard tools for the fluctuations analysis, both experimentally as well as numeri¬ 
cally, are the nearest-neighbor spacing distribution P{s) and the number-variance S^(r), the 
measures for the short and long-range spectral correlations, respectively ElE]. These can be 
dehned as follows: P{s) (same as P{l;s,e) defined in section II. A) is the probability of two 
nearest neighbor eigenvalues to occur at a distance s, measured in the units of local mean 
level spacing D, and, S^(r) is the variance in the number of levels in an interval of length 
rD. Figure 7 compares P{s) for Gaussian and bimodal disorders for various dimensions for 
hxed size N = 10^ and for two energy-ranges i.e the edge and bulk of the spectra; it clearly 
shows the sensitivity of P{s) to both dimensionality as well as the energy-range but almost 
no effect of the nature of disorder i.e Gaussian or bimodal. (note, as hgures 7(a,d) indicate, 
the nature of disorder seems to have some influence only for <7=1). For lower dimension 
or in the edge region, the statistics is shifted more towards Poisson, an indicator of the 
increasing eigenfunction localization. For <7=1, the P{s) behavior seemingly approaches 
singularity at s = 0, implying a high degree of quasidegeneracy among eigenvalues. As, for 
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d = 1, the system is nearly integrable, the singular behavior is in qualitative accordance with 
Shnirelman theorem which states that for a classically near integrable system at least each 
second level spacing in the corresponding system become exponentially small [13]. (More 
formally, the theorem states that the spectrum is asympotically multiple i.e for each m > 0 
there exists Cm > 0 such that min(efc — Ck-i, Ck+i — eu) < [T3]). 

To analyze size-dependence of the short range correlations, we consider P{s) behavior 
for four different sizes of d = oo (hgure 8 ) and d = 3 cases (hgure 9). Although sensitive 
for small values, the distribution shows a tendency to become size-independent for large N- 
values. As the hgures indicates, the limiting distribution for both dimensions is independent 
of the nature as well as strength of disorder; this is in accord with our theoretically claimed 
Gaussian-bimodal analogy discussed in section 4 of [T| (between cases II and IV for d = cx) 
and cases III, V for hnite d). 

As mentioned above, the long-range fluctuations around the spectral density can be 
analyzed by S^(r). Figures 10, 11 show S^(r) behavior of the CCE spectra for both bulk 
and edge regions, of d = oo as well as d = 3 cases, respectively; the hgures reconhrm the 
behavior indicated by P(s)-numerics discussed above i.e independence from the nature or 
strength of disorder but sensitivity to basis-connectivity and/ or energy-range. The S^(r) 
numerics also indicates large V-approach to size-invariant behavior for both d = oo, 3 cases 
but the convergence is slower in the edge region and for large r-values. (Note, the curves 
for N = 512 in hgure 10 and hgure 11 seem to display a non-increasing behavior but as 
mentioned in section II.C, the long range measures contain spurious huctuations for small 
system sizes and true huctuations emerge only in inhnite size limit | 6 |. The numerics for long 
range measures therefore can not be considered reliable by itself however it can be used as 
a support for the theoretical results or for other measures (e.g. P{s)). As clear from hgures 
10,11, the curves for large N not only show smaller deviation from each other (e.g. N = 2197 
and N = 4913) (indicating approach to size-independence as predicted theoretically by us) 
but their behavior is also consistent with corresponding P{s) behavior. Also note, due to 
hnite size ehects, the statistical error in S^(r)-analysis increases with large r). 
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III. CCE-BE ANALOGY 


The theoretical analysis discussed in section IV of [T] indicates that a N x N CCE of real- 
symmetric matrices H is analogous to a section of the BE i.e the sub-ensemble with uniform 
eigenvector and zero eigenvalue. As eq.(72) of P indicates, the local spectral statistics of a 
N X N CCE approaches, in large N limit, the BE-statistics. The ensemble density for the 
BE analog can be given as 


p{H) oc exp 


R 

2 


N 




N 

r,(l + cAf) 


( 6 ) 


with c = 1 for the analog of a CCE with d = oo] (note the theoretical analysis in P could 
not exactly predict c-value for the BE analog of CCE, finite d case). Eq.([^ is also known as 
Rosenzweig-Porter (RP) (or Porter-Rosezweig) ensemble [15]. Referring it as a BE, however, 
has some advantage: it gives an incentive to seek connections among BEs and CCEs under 
global symmetries other than time-reversal. 

As expected due to non-stationarity, the local fluctuations of a BE analog vary along the 
spectrum axis. Due to relatively smaller level-spacing, the statistics in the bulk of the BE 
spectrum is very close to a GOE but the edge-statistics deviates signihcantly from the latter 
ng. This indicates that the GOE-CCE analogy is expected to be valid only for a very small 
spectral range AE around e = 0 P which becomes narrower as the size N increases. The 
theoretical study in P however claims that the 2-point level-density correlation R 2 {Ti,r 2 ) 
{the probability of finding two levels at a distance |ri — r 2 |) for a GGE is analogous to that of 
a GOE for all energy ranges. To clarify this discrepancy and hnd the correct GGE analog, 
we numerically compare the spectral GGE statistics with that of a BE and a GOE, both 
in the bulk as well as edge. Figure 12 displays the comparative behavior for four spectral 
measures, namely, Ri{r), R 2 {ri,r 2 ) as well as P{s) and S^(r) for GGGE, d = oo and BE, 
each consisting of 5000 matrices of size N = 10^. As can be seen from figure 12(a), i?i(e) 
for the d = oo GGGE case deviates a little from its BE analog but their local fluctuations 
measured in terms of P{s), S^(r) and R 2 {f') are analogous. (Although R 2 {r) provides the 
same information as given by P{s) and S^(r)), it is considered here to substantiate our 
analytical claims against the one in P). For the clarity of presentation, the d = oo bimodal 
case is not shown in figure 12 but its analogy with corresponding GGGE case (figures 7-9) 
implies the same with BE. The above comparison is also repeated for GGGE, d = 3 case; 
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as confirmed by figure 13, its BE analog, obtained by a numerical search, corresponds to 
fi ^ 200N. 

In general, the BEs represent the non-equilibrium states of transition, driven by a single 
parameter, between two universality classes of random matrix ensembles. But the BE analog 
of a CCE with inhnite connectivity is independent of any parameters except the system-size 
N. In asymptotic limit iV —)■ oo, therefore, the statistical fluctuations of a CCE with d = oo 
are free of all parameters and represent a new universality class (different from ten standard 
universality classes). (As seen in hgures 8 and 10, both P{s) and E^(r) for CCE, d = oo 
are size-independent; note, in hgure 10, an apparent size-dependence of S^(r) for small N 
is due to spurious fluctuations introduced by hnite size effects). 

IV. CONCLUSION 

In this paper, we have presented a detailed numerical analysis of the spectral statistics 
of the column constrained ensembles (real-symmetric case with independent off-diagonals 
and column constant a = 0 for all columns) under various system conditions. Our study 
clearly indicates a non-ergodic as well as non-stationary behavior of the CCE-spectrum 
even for inhnite connectivity; this is in contrast to the ergodic, stationary nature of GOE 
(which is similar to the CCGE in terms of matrix structure except for column constraints). 
This implies an energy dependence of the spectral huctuations which is conhrmed by our 
numerics. We also hnd that the strength or nature of disorder has no effect on the local 
spectral-huctuations. Although the average behavior of both, level density as well as inverse 
participation ratio, is sensitive to the type of disorder (for hnite basis-connectivity), its de¬ 
pendence on the disorder-strength can be scaled out. (Note here only those types of disorder 
are implied which do not change the global constraint class of the ensemble). Further both 
the average behavior as well as the huctuations are sensitive to other system conditions e.g. 
dimensionality/ basis-connectivity and energy-range. As connectivity in the basis increases, 
the average behavior also become insensitive to nature of disorder. 

In large-size limit, the spectral huctuations in CCEs become size-independent too. As 
predicted in [T] for CCE case with inhnite connectivity, our numerics conhrms the analogy 
of its huctuations to those of a special type of critical BE, intermediate between Poisson 
and COE and free of all parameters. But, for a CCE with hnite connectivity, the BE analog 
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depends on a single parameter. This indicates a cross-over of CCE-statistics from Poisson 
to this special critical BE with basis-connectivity as the transition parameter. 

The present study still leaves many questions unanswered. For example, the real- 
symmetric ensembles considered here are applicable to time-reversal systems with column 
constraints. It would be interesting to numerically analyze the complex Hermitian ensem¬ 
bles with column/row constraints; these ensembles are appropriate models, for example, for 
Goldstone systems without time-reversal symmetry). The presence of correlations among 
off-diagonals (applicable to interacting systems), varying column constants for each column 
or an absence of Hermiticity (relevant for systems e.g Google matrix) are some other impor¬ 
tant effects which are yet to be probed. Another important open question is the effect of 
column constraints on the eignvector fluctuations. It would also be interesting to explore the 
exact dependence of the GGE-statistics on the relative values of the ensemble parameters 
and the column constraints e.g ratio a /^/7 in case III of [1]. We expect to answer some of 
these questions in near future. 
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FIG. 1: Ergodic behavior of level density Psm{&)'- The figure compares the spectral averaged 
level-density Psm{e) for a single matrix to an ensemble average i?i(e) as well as {psm{&)) for four 
ensembles: (a) CCGE (eq.Q), (b) CCGE (eq.Q, (c) BE (eq.Q), (d) GOE. In each case, we 
consider an ensemble of 5000 matrices of size N = 1000. The numerics shows that psm for CCGE 
and BE fluctuates from one matrix to the other which manifests in a clear, although small, deviation 
of {psra) from Ri] the deviation however seems to reduce with increasing N. As shown in the part 
(d), the ergodic nature of GOE is clearly visible from an exact analogy of i?i(e) with {psm) and 
also with Psm for a single matrix (to emphasize the latter, two matrix cases are shown in figure 
1(d)). 
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FIG. 2: Non-ergodicity of local fluctuations: The figure compares the spectral averaged 
behavior of with ensemble average for many k values and for spacings taken from bulk near 
e = 0; here measures the width of the eigenvalue spacing distribution P{k,x,e), to find k^^ 
nearest neighbor levels at a distance x in a spectral range around e. The spectral and ensemble 
averages are obtained by using eq.(36) and eq.(39) of [6] respectively. The results are shown for 
four different ensembles: (a) CCGE, d = oo, eq.Q, (b) GGGE, d = 3, eq.Q, (c) BE, eq.([^, (d) 
GOE. The deviation of ensemble average from the spectral one in first three case is an indicator of 
non-ergodicity of local fluctuations which however seems to be quite weak. The non-ergodicity is 
confirmed by the lack of stationarity of shown in hgure 3. Although a small deviation is seen in 
the case of a GOE too but it is due to spuriouslffuctuations and is a hnite size effect; (as indicated 
by the stationarity displayed in figure 3(d)). 
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FIG. 3: Stationarity of local fluctuations: The figure displays the behavior of ensemble- 
averaged CTfc, the width of the eigenvalue spacing distribution p(k,s,e) for many k values and for 
20 level-spacings taken from the levels of five different spectral-regimes: (i) AEi = 98 — 118, (ii) 
AE 2 = 196 - 216, (hi) AE 3 = 294 - 314, (iv) A.F 4 = 392 - 412, (v) AE 5 = 490 - 510, Again 
four ensembles are considered: (a) CCGE, d = 00 , eq. 0. (b) CCGE, d = 3, eq. (©, (c) BE, eq.@, 
(d) GOE. Note the deviation of cr^ for a fixed k and for different spectral-ranges is significant in 
parts(a,b,c) while it is almost negligible in part(d); this is an indicator of the non-stationarity for 
the first three cases. 
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FIG. 4: Disorder-dependence of level density Psm{^)'- The figure illustrates the behavior 
of Psm{.&) of a column-constrained ensemble for a fixed size N = 1000 but for different disorder- 
strengths 7 = 2 ^ for Gaussian disorder and different basis-connectivity (to consider dimensionality 
effects): (a) Gaussian case d = 1, eq.Q (b) Gaussian case d = 3, eq.Q (c) Gaussian case, d = oo, 
eq.Q. Note, in parts (a,b,c), psm is scaled by ^ which results in an analogous form for 

different cr; this scaling is same as used in eqs.(34,54) of [T] to remove 7 -dependence. To understand 
the dependence on the disorder-type, we also consider bimodal case for a fixed disorder strength 
and three types of basis-connectivity (for N = 1000): (d) bimodal case, d = 1, eq.Q (e) bimodal 
case, d = 3, eq. 0 (f) bimodal case IV, d = 00 , eq. 0 . 
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FIG. 5: Size-dependence of level density psm{e)'- figures (a)-(c) describes the level-density 
behavior of a column-constrained ensemble for a hxed disorder-strength = 1 ( 7 = ( 2 cj^)“^)) 
and for different sizes N-. (a) d = 1, Gaussian case eq.Q, bimodal case eq.Q, (b) d = 3, Gaussian 
case eq.Q and bimodal case eq.Q, (c) d = 00 , Gaussian case eq. Q and bimodal case eq.Q. 
To avoid cluttering, the notations in (a)-(c) are changed; the symbols Gi,G 2 ,G 3 ,G 4 refer to 
Gaussian cases with N = 512,1000,2197,4913 respectively and symbols Bi, B2, B3, B4 refer to 
Bimodal cases with N = 512,1000,2197,4913 respectively. Note in parts(a,b) psm is rescaled: 
Psm —^ Psm/N but in part (c), the analogy of different N-cases results only after following rescaling: 
e —>■ e/y/N, Psm —^ Psm/VN. As shown later in figure 12(a), psm for d = 00 case behaves as a 
semi-circle in the bulk: psm{e) = — e^. The parts (d)-(f) depict the density of the low- 

energy excitations for different dimensions and two disorder-types, along with fits /i for Gaussian 
and /2 for bimodal case: (d) d = 1, /i = 0.003 , f 2 = O.Olx*’'^, (e) d = 3, /i = 7.5 x 10“^ 

/2 = 0.004 x^, (f) Infinite range case, /i = /2 =201.6 x 10“^x^/^. The fits in hgures (e)-(f) seem to 
agree with available results for Goldstone modes in spin glasses. 
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FIG. 6: System-dependence of inverse participation ratio l 2 '. The figure describes the 


/ 2 -behavior of a column-constrained ensemble {N = 1000) for many disorder strengths (with 7 = 


{2a 


2\-l^ 


(only for Gaussian disorder), basis-connectivity and for two different type of disorders. 


The insets in each case show the behavior for small excitation-energy u along with the fits by solid 
lines: (a) Gaussian case eq.Q, d = 1, fit = 0.479 — 0.0028x*^'^, (b) Gaussian case eq. (§, d = 3, 
fit = 0.485 —0.028x^/^ — O.Olx^, (c) Gaussian case eq.Q, d = 00 , fit = 0.61 —0.025x-|-5 x 10“® x^, 
(d) bimodal case eq.Q, d = 1, fit = 0.05 -|- 0.32a;°'^, (e) bimodal case eq.Q, d = 3, fit = 
0.355 —O.lx^/^-1-0.3x^ (f) bimodal case eq.([^, d = 00 , fit = 0.61 — 0.028x4-2 x 10“^ x^. The plots 
in (a, b, d, e) are rescaled by the off-diagonal variance The infinite-range case for Gaussian is 
considered for a single disorder-strength only; (as clear from eq.Q, the off-diagonal variance for 
this case can be scaled out). 
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FIG. 7: Dimensionality-dependence of nearest-neighbor spacing distribution: The figure 
compares P(s)-behavior for two different types of disorders, and three types of basis-connectivity 
of the Column constrained ensemble in two different energy ranges. The size is kept fixed at 
N = 1000; note, N = L'^ for cases with d = 1,3. The Gaussian and bimodal cases for d = 3 
corresponds to eq. (§ and eq. @ respectively . For d = oo, eq. 0 and eq. 0 give the Gaussian 
and bimodal cases respectively. The parts (a)-(c) describe the statistics in the edge and (d)-(f) in 
the bulk of the spectra: (a) edge, d = 1, (b) edge, d = 3, (c) edge, d = oo, (d) bulk, d = 1, (e) 
bulk, d = 3 (f) bulk, d = oo. The figures suggest an increase of level-repulsion with increasing 
basis-connectivity, with rate of change slower in the edge than bulk. The difference in the edge 
and bulk statistics for each dimensionality also indicates a lack of stationarity in the spectrum. 


22 




























s 



s 



s 



s 


FIG. 8: Size-dependence of P{s) for colnmn constrained ensemble with d = oo: The 

figure describes the P(s)-behavior for the Gaussian case (eq.Q and bimodal case (eq.Q) for four 
system sizes N. As seen in the hgure, P{S) approaches to an invariant form as the system size N 
increases. The limit is sensitive to the energy-range (e.g. edge vs bulk) but almost independent of 
the nature of disorder (e.g. Gaussian vs bimodal) :(a) edge, Gaussian, (b) edge, bimodal, (c) bulk, 
Gaussian, (d) bulk, bimodal. Note the statistics in the bulk is close to a GOE as expected for a BE 
with jjL = N (only 2% levels chosen from the center of the bulk in each case). In the edge regime, 
intermediate behavior of P{s) to Poisson and GOE limits is again in agreement with theoretically 
expectation (as ~ 1 > A^dge (see section V.B of [I])). 
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FIG. 9: Size-dependence of P{s) for d = 3 case: The figure shows the P(s)-behavior for 
Gaussian case (eq. Q and bimodal case V (eq. Q for four sizes N for the case d = 3 where N = 
and periodic boundary conditions are imposed at length L. Here again the approach of P{S) to 
an invariant form with increasing N is indicated. Similar to d = oo case, here also the approach 
is sensitive to the energy-range (e.g. edge vs bulk) but is almost independent of the nature of 
disorder (e.g. Gaussian vs bimodal) : (a) edge, Gaussian, (b) edge, bimodal, (c) bulk, Gaussian, 
(d) bulk, bimodal. As compared to d = oo case (fig 8), the edge-statistics here is shifted more 
towards Poisson regime which is expected due to an increased localization of the eigenfunctions, 
originating in confined hopping. 
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FIG. 10: Size-dependence of the number-variance S 2 (r) for case d = oo: figure shows the 
variance of number of levels in a distance of r mean level spacings for Gaussian case (eq. 0), and 
bimodal case (eq.0): (a) Gaussian, edge, (b) Bimodal, edge, (c) Gaussian, bulk, (d) Bimodal, 
bulk. To detect the A^-sensitivity in —>• oo limit, each case is considered for four system sizes. 
As visible in parts 10(c,d), the statistics in the bulk is again near a GOE and reconfirms P{s) 
behavior shown in figure 8(c,d). Note, the approach to size-invariance for S^(r) in the edge region 
is slower for large r but is visible from the behavior of curves for N = 2197,4913 in parts 10(a,b). 
The statistics now is intermediate between Poisson and GOE limits; this is in conformity with 
corresponding P{s) behavior shown in figure 8(a,b). 
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FIG. 11: Size-dependence of number-variance S 2 (r) for d = 3 case: The Figure shows the 
variance of number of levels in a distance of r mean level spacings for the Gaussian case (eq.Q), 
and, the bimodal case (eq.Q for four sizes N: (a) Gaussian, edge, (b) Bimodal, edge, (c) Gaussian, 
bulk, (d) Bimodal, bulk. Here again, the approach to size-invariance in the edge region is slower 
for large r but is visible from the behavior of large A^-curves; similar to P{s) behavior for d = 3 
(figure 9), the statistics here is intermediate between Poisson and GOE limits. 
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FIG. 12: Comparison of column constrained ensemble, d = oo with BE: Here we compare 
CCGE given by eq.Q along with its BE analog (given by eq.Q with // = N). Each ensemble is 
considered for a fixed size N = 1000 but the results shown are independent of the system size in 
the study: (a) Ri{e){ rescaled by N), (b) i? 2 (e), edge, (c) R 2 {e), bulk (d) P{s), edge, (e) P{s), 
bulk, (f) S^(r), edge, (g) S^(r), bulk. The part (a) clearly indicates the deviation, although small, 
of Ri{e) for GCGE from that of BE; as expected, the semi-circle fit= ^ 3 ^jy V agrees well 
with the bulk of BE but deviates in the tail regime. As seen in the parts (b,c), for CCGE 

deviates from that of GOE in the edge but agrees well in the bulk. The parts (d)-(g) confirm the 
CCE-BE analogy of the local fluctuations for (?2ussian disorder, both in the bulk as well as in the 
edge-region. The same analogy is valid for CCE case with bimodal disorder too which, although 
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FIG. 13: Comparison of column constrained ensemble, d = 3 case with BE: Here we 
compare CCGE given by eq.([^ with the BE, given by eq.Q with c = 200. Each ensemble is 
considered for a fixed size N = 1000 and four four fluctuation measures: (a) Ri{e){ rescaled by 
N), (b) i? 2 (e), edge, (c) R 2 ie), bulk (d) P{s), edge, (e) P(s), bulk, (f) edge, (g) ^^(r), bulk. 

The deviation of i?i(e) for the GGGE from that of BE is now clearly visible from the part(a). But, 
as the parts(b, c, d, e) indicate, R 2 and P{s) in both bulk and edge agree well for the two cases. A 
small deviation seen in the parts (f,g) for S 2 (r) may be a finite size-effect. Another possibility is 
that the BE with c = 200 is although close but^ not an exact analog of the CCGE case considere 
here. Note our theoretical analysis given in section IV.B of [T] does not exactly predict the BE 








































